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Abstract 

The long-time behaviour of solutions of systems of conservation laws has been extensively studied. In particular, Liu and 
Zeng O have given a detailed exposition of the leading order asymptotics of solutions close to a constant background state. 
In this paper, we extend the analysis of [4] by examining higher order terms in the asymptotics in the framework of the 
so-called two dimensional p-system, though we believe that our methods and results also apply to more general systems. 
We give a constructive procedure for obtaining these terms, and we show that their structure is determined by the interplay 
of the parabolic and hyperbolic parts of the problem. In particular, we prove that the corresponding solutions develop long 
tails that precede the characteristics. 



1 Introduction 



In this paper, we consider the long-time behavior of solutions of systems of viscous conservation laws. This topic has been 
extensively studied. In particular, for the case of solutions close to a constant background state, [4| contains a detailed 
exposition of the leading order long-time behavior of such solutions. More precisely, it is shown in |4| that the leading 
order asymptotics are given as a sum of contributions moving with the characteristic speeds of the undamped system of 
conservation laws and that each contribution evolves as either a Gaussian solution of the heat equation or as a self-similar 
solution of the viscous Burger's equation. Thus with the exception of the translation along characteristics, these leading 
order terms reflect primarily the dissipative aspects of the problem. 

In this paper, in an effort to better understand the interplay between the hyperbolic and parabolic aspects of the problem, 
we examine higher order terms in the asymptotics. We work with a specific two-dimensional system of equations - the 
p-system, but we believe that its behavior is prototypical. In particular, we think that our methods and results would extend 
to more complicated systems such as the 'full gas dynamics' and the equations of Magneto-Hydro-Dynamics (MHD) as 
considered in @). 

The specific set of equations we consider is the following: 

d t a = c 1 d x b, a(x,0) = a (x) , 

d t b = c 2 d x a + d x g(a, b) + a (d x b + d x (f(a,b)d x b)) , b(x, Q) = b (x) . 

We will make precise the assumptions on the nonlinear terms / and g below, but in order to describe our results informally, 
we basically assume that \g(a, b)\ — 0((\a\ + \b\) 2 ) and \ f(a,b)\ ~ C((|a| + We also note that without loss of 
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generality, we can set c\ = Ox = 1 and a = 2 in ( 11.11 ). which can be achieved by appropriate scalings of space, time and the 
dependent variables, and possible redefinition of the functions / and g. 

Physically, fll.lt is a model for compressible, constant entropy flow, where a represents the volume fraction (i.e. the 
reciprocal of the density) and b is the fluid velocity. The first of the two equations in fll.ll i is the consistency relation 
between these two physical quantities. In particular, it would not be physically reasonable to include a dissipative term in this 
equation, whereas such a term arises naturally in the second equation which is essentially Newton's law, in which internal 
frictional forces are often present. As a consequence of the form of the dissipation the damping here is not 'diagonalizable' 
in the terminology of (4) . 

Next, we note that with the scaling ci = C2 = 1 and a = 2 in d 1 - lb - the characteristic speeds are ±1. Then, following 
Liu and Zeng O, we introduce new dependent variables u and v which translate with those characteristic speeds ±1, 
respectively. If the initial conditions ao and bo in ( 11. lb decay sufficiently fast as |x| — > oo, Liu and Zeng showed that in the 
translating frame of reference, u(x,t) — ^ +t 9o( ^fjp ) + C((T + and similarly for v, where go is a self-similar 

solution of either the heat equation, or of Burger's equation, depending on the detailed form of the nonlinear terms. In this 
paper we derive similar expressions for the higher order terms in the asymptotics through a constructive procedure that can 
be carried out to arbitrary order. 

More precisely, we show that for any N > 1, there exist (universal) functions {<7„ }„Li and constants {d^}^ =1 determined 
by the initial conditions, such that 

u(x,t) = —=g+(-£-) + £ x j + ( . (1.2) 

We give explicit expressions for the functions below, but focusing for the moment on the case N = 1 and the variable 
u, we have 

u(x, t) = —^g+t-JL-) + L_^ fl +(* ) + o( ^-r) , 

K ' VTTt Joy ^ I + i ' (i+t)i iyiv vr+* ; V(i+t)i;' 

where the functions g^(z) and gf(z) are solutions of the following ordinary differential equations: 

dldt^) + \zd z g+(z) + i 9o +(z) + c + d z (g+(z) 2 ) - (1.3) 

d 2 z gt(z) + \zd z gf(z) + - A g+{z) + 2c+d z (g+ (z)g+ (z)) = . (1.4) 

Here c + is a constant that depends on the Hessian matrix of g(a, b) at a = b = and that will be specified in the course 
of our analysis. We will prove that while all solutions of ( 11. 3b have Gaussian decay as \x\ — > oo, general solutions of the 
linear equation ( 11.41 ) are linear combinations of two functions g^ ± (z), where g^ ± (z) decays like a Gaussian as z — > =Foo 

but only like |z| - ^ as z — > ±oo. The graphs of the functions <?j^(z) anc ^ 5i" ( z ) ^ presented in FigureQ] 

Thus, the higher order terms in the asymptotics develop long tails. These tails are a manifestation of the hyperbolic part of 
the problem (or perhaps more precisely of the interplay between the parabolic and hyperbolic parts). Were we to consider 
just the asymptotic behavior of the viscous Burger's equation which gives the leading order behavior of the solutions, we 
would find that if the initial data is well localized, the higher order terms in the long-time asymptotics decay rapidly in space 
and have temporal decay rates given by half-integers. 

Another somewhat surprising aspect of our analysis is that the tails actually precede the characteristics. 

We also note one additional fact about the expansion in (11.2b . Prior research J2] [7] has shown that for both parabolic 
equations and damped wave equations the eigenfunctions of the operator 

Cu(z) — d\u + -zd z it 
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play an important role for the asymptotics. In particular, on appropriate function spaces this operator has a sequence of 
isolated eigenvalues whose associated eigenfunctions can be used to construct an expansion for the long-time asymptotics. 
In this connection we prove that the functions are closely approximated by eigenfunctions of C with eigenvalues X n = 
— j + 2~ ; more precisely, the functions are eigenfunctions of a compact perturbation of C, see e.g. dl.41 >. However, 
so far we have not succeeded in finding a function space which both contains these eigenfunctions (the functions decay 
slowly as z — > ±00) and in which the corresponding eigenvalues are isolated points in the spectrum. We plan to investigate 
this point further in future research. 

Before moving to a precise statement of our results we note that our approach makes no use of Kawashima's energy estimates 
for hyperbolic-parabolic conservation laws [3 1. Instead we prove existence by directly studying the integral form of (ll.lt . 

We now state our results on the Cauchy problem (II . lb - We begin by stating the precise assumptions we make on the 
nonlinearities / and g in (II. lb . 



Definition 1 The maps /, g : R 2 — > R are admissible nonlinearities for ( 17.71 ) if there is a quadratic map go : R 2 — > R and 
a constant C such that for all |z|, |zj| and | Z2 1 small enough, 

| ff (z)| < C|z| 2 , | 5 ( Zl ) - g(z 2 )\ < C|z x - zalflzxl + |z 2 |) , 

|A ff (z)| < C|z| 3 , |A 5 ( Zl ) - A 3 (z 2 )| < C\l x ~ Z2|(N + |z 2 |) 2 , 
|/(z)|<C|z| and |/( Zl )-/(z 2 )| <C| Zl -z 2 | , 

where Ag(z) = g(z) ~ go(z). 



The main result of this paper can be formulated as follows: 



Theorem 2 Fix N > 0. There exists eo > sufficiently small such that if 
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(i) |ao|H!(R) + holi 1 ^) < eo an d |6o|h 2 (k) + l^ol^ 1 ^) < eo 

(ii) \x 2 a \ L 2 m + |a; 2 6o|^ (R) < oo, 



then ( 17.71 ) has a unique (mild) solution with initial conditions ao and bo- Moreover, there exist functions {<?„}^=o (indepen- 
dent of initial conditions) and constants Cjy, {d„}„ = i determined by the initial conditions such that if we define 



u(x, t) = a(x — t,t) + b(x — t, t) and v(x, t) — a(x +t,t) — b(x + t, t) 



then 



u(x, t) 
v(x, t) 



1 N 1 

^ + (^=) + E 7; — o=z^<fln (ttr) + R u^ *) 



„=1 (l + O^^r 



where the remainders and 7?^ v satisfy the estimates 



N 



(1-5) 



S U V (l + t)*-^T*\\R? UtV} (;t)\\ L z < C N 



t>0 



sMi + t)*-^\\d x R^ v} (-,t)\\L 2 ( m ^ c * 



(1.6) 



Furthermore, for n > 1, the functions g^ satisfy g^(z) ~ |z| 



-1+2- 



fliZ-> ±00. 



There is a slight incongruity in this result in that the norm in which we estimate the remainder term is weaker than that we 
use on the initial data; namely, we do not give estimates for the remainder in H 2 (R), or in the localization norms L 1 (K) 
and the weighted L 2 (R)-norm (on that aspect of the problem, see Remark [3] below). Theorem factually holds for slightly 
more general initial conditions than those satisfying (i)-(ii). Furthermore, we will prove that the estimates (11.61 ) hold for all 
initial conditions (ao, bo) in a subset T>2 C Hi x H2 that is positively invariant under the flow of ( ll.lt . However, since the 
topology used to define the subset T>2 is somewhat non-standard, we have chosen to state the result initially in this slightly 
weaker, but hopefully more comprehensible, form to keep the introduction as simple as possible. 



Remark 3 It is interesting to note (see Proposition\7\below) that ||cc 2 a(-, ^H^ 2 ^ + ll a;2 &('; ^) Hi 2 (r) is finite for all finite 
t > 0, but that the terms with n > 1 in the asymptotic expansion do not satisfy this property due to the long tails of the 
functions g^. 



Remark 4 As the asymmetry in the degree of x derivatives in ( I7.7D suggests, we require more spatial regularity from the 
second component (the b variable) than from the first (the a variable). It is then natural to expect that R 1 ^ or R 1 ^ are not 
necessarily in H 2 , but that only their difference is. 



We conclude this section with a few remarks. Define u± (x, f) = a(x, t) ± b(x, t). Then the asymptotics of the solutions of 
dl.lt in the variables u± are the same as those of the two dimensional (generalized) Burger's equation 



dtu+ = d x u + + d x u + + d x (c+u + — c_u_) 
9 f u_ = <9 2 u_ — d x ii_ + d x (c_u_ — c + u 2 + ) , 



(1.7) 



where the constants c± are determined by the Hessian of g(a, b)ata — b — through 



c±=±-(l,±l) 



dig d a d b g 
d a d b g dig 



=b=Q 



1 

±1 
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We will see that the hyperbolic effects manifest themselves through the 'source' terms — C-V?_, respectively c + u\ in the 
first, respectively second equation in ( 11.7b . In particular, none of the terms with n > 1 would be present in the asymptotic 
expansion if those terms were absent. 

Finally, note that we have chosen to state Theorem [2] for finite N. As it turns out, the sums appearing in jl.5\ converge 
in the limit as TV — > oo, in which case the estimates ( 11.61 ) hold with time weights replaced by (1 + t)* ln(2 + t)~ x and 
(1 + t) 3 ln(2 + i) -1 . The proof can easily be done with the techniques used in this paper and is left to the reader. 

The remainder of the paper is organized as follows: In Section [2] we discuss the well-posedness of the Cauchy problem 
( II. U in an appropriately defined topology. In Section [3] we explain our strategy for proving our main result, Theorem [2] 
on the long time asymptotics of solutions of ( II. lb . Namely, we decompose that proof into a series of simpler sub-problems 
which are then tackled in subsequent sections: in Sections H] and [5] we investigate properties of solutions of Burger's type 
equations, respectively of inhomogeneous heat equations, as they occur naturally in the asymptotic analysis. In Section [6] 
we collect some estimates that are used in the proof of the well-posedness of ( 11.11 ). Finally, in Section [TJ we specify the 
sense in which the semigroup of the linearization of ( II. U is close to heat kernels translating along the characteristics, and 
we give estimates on the remainder terms occurring in Theorem[2] 



2 Cauchy problem 



To motivate our technical treatment of the problem and in particular our choice of function spaces, we first note that upon 
taking the Fourier transform of the linearization of dl.ll i, it follows that 



v&/ \bj ~\ik -2k t 
We then find that the (Fourier transform of) the semigroup associated with ( 12.11 ) is 

L< _ -k 2 t( cos(fctA) + 2j sin(HA) ± sin(HA) \ 

V i sin(£:iA) cos(fciA) - | sin(fctA) J ' y ' 

where A = \/l — k 2 . The most important fact about the semigroup e is that it is close to c^ ot , the semigroup associated 
with the problem 



(), 



Formally, e^°* can be obtained by setting A = 1 in e^* and by conjugating with the matrix 

These two operations correspond to a long wavelength expansion and a change of dependent variables to quantities that 
move along the characteristics. More precisely, we will prove that e^ 4 satisfies the intertwining property 

Se Lt w e Lot S , 

where the symbol rj means that the action of these two operators is the same in the large scale - long time limit; see Lemma 
[T9lat the beginning of Section|7]for details. 

Furthermore, e^* satisfies parabolic-like estimates 



e Lt | ^Ce-^-^n ^+^Y (2.5) 
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(2.6) 



uniformly in t > and ieR. 

Hence, to summarize, er"* behaves like a superposition of heat kernels translating along the characteristics of the underlying 
hyperbolic problem. In view of the above observations as well as of classical techniques for parabolic PDE's, see e.g. J5][T1, 
we will consider ( II. U in the following (somewhat non-standard) topology (cf also @): 

Definition 5 We define Bo, resp. B, as the closure q/C^°(R, R 2 ), resp. CJ°(R x [0, oo), R 2 ), under the norm \ ■ \, resp. 
| • ||, where for zq = (do, bo) : R — ► R 2 and z = (a, b) : R x [0, oo) — * R 2 , we define 

|z | = Hzolloo + II z 1| 2 + ||Dz || 2 + ||D 2 6 || 2 , ||z|| = Hzll^o + ||z|| 2i i + ||Dz|| 2) | + ||D 2 6|| 2!r . 

Here (Da)(x, t) = d x a(x, t), a(k, t) is the Fourier transform of a(x, t), 

\\f\\ M = SU P (1 + ty\\f(;t)\\ P , ll/IU* = SU P ^±^-\\f(;t)\\ P 
t>0 t>0 ln \ z + t) 

and || • ||p is the standard L p (R) norm. 

Before turning to the Cauchy problem with initial data in Bo we collect a few comments on our choice of function spaces. 

Consider first the requirements on the initial conditions in dl.lt . While the use of H 1 space is quite natural in this context, 
we choose to replace the L 1 norm by the (weaker) control of the L°° norm in Fourier space. This has the great advantage 
that all estimates can then be done in Fourier space, where the semigroup e^* has the simple, explicit, form M.2\ . 

In turn, our choice of ^-exponents in the norm || • || is motivated by the fact that these are the highest possible exponents for 
which the || • ||-norm of the leading order asymptotic term ^\ +t 9o{ 7/f!pf ) * s bounded. Note also that for the linear evolution 
(12. U . we have 

||e Lt z || < C|z | , (2.7) 

since j(k,t) = e~ min (' c2 ' 1 ^uo(fc) satisfies 

||D n j(- > *)||a<a(e- t ||D n «o||a + miu(«-*-* ||«o||oo, |p n «o|| 2 )) 

for all n — 0,1, . . .. 

Finally, we note that for admissible nonlinearities in the sense of DefinitionQ] the map h(a, b) = f(a,b)d x b+g(a,b) = h(z) 
satisfies 

\\h(z)\\ lti + ||ft(z)|| 2 ,3 + \\Dh(z)\\ 2<i < C\\z\\ 2 , (2.8) 
\\h{zx) - h{z 2 )\\ lt t + - /i(z 2 )|| 2) | < C||z x - z 2 ||(|| Zl || + ||z 2 ||) , (2.9) 

\\D{h{z x ) - ft(z 2 ))|| 2i 5 < C||z x - z 2 ||(|| Zl || + ||z 2 ||) . (2.10) 

We are now fully equipped to study the Cauchy problem ( 11.1b in B: 

Theorem 6 For all zq € Bo with \zq\ = \(ao, bo)\ < eo small enough, the Cauchy problem ( li.iD is (locally) well posed in 
B if the nonlinearities are admissible in the sense of Definition^ In particular, the solution satisfies ||z|| < ceo for some 
c > 1 and is unique among functions in B satisfying this bound. 
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Proof. Upon taking the Fourier transform of dl.lt . we get 

which gives the following representation for the solution 

*«> S Q - eL '0 (.))) S +Af W(0 ■ (2-12) 

We will prove below that for all Zj € 0, i = 1, 2, we have 

||W[z]|| < C||z|| 2 and ll^zi]-^]!! <C||zi-z 2 ||(||zi|| + ||z 2 ||) (2.13) 

for some constant C. The proof of Theorem [6] then follows from the fact that for all zq G Bq with |zq| < eo small enough 
and c > 1, the r.h.s. of ( 12.121 ) defines a contraction map from some (small) ball of radius ceo m B onto itself. 

The general rule for proving the various estimates involved in ( I2.131 l is to split the integration interval into two parts, with 
s eli e [0, |] and s e I2 e [|, t]. In 2i, we place as many derivatives (or equivalently, factors of k) as possible on the 

semigroup e (* _s \ while on Z 2 , (most of) these derivatives need to act on h, since the integral would otherwise be divergent 
at s = t. 

Additional difficulties arise from the fact that e^' has very little smoothing properties (slow or no decay in k as \k\ — > 00), 
so that in some cases we need to consider separately the large-fc part and the small-fc part of the L 2 norm, say. This is done 
through the use of P, defined as the Fourier multiplier with the characteristic function on [—1, 1]. 

We decompose the proof of ||A/"[z] |j < C||z|| 2 into that of 

\\AT[z]\\ < HATMIUo + ll-^N|| 2 .i + ||MW[z]|| a 3 + 11(1 - P)DAA[z]|| 2 3 

+ ||(1 - P)D 2 AA[z] 2 || 2ir + ||(1 - Q)PD 2 AA[z] 2 || 2jr + ||QPD 2 AA[z] 2 || 2 , r 
<C||z|| 2 , (2.14) 

where Q is the characteristic function for t > 1 and jV[z] 2 denotes the second component of A/"[z]. 
We now consider ||PDjV[z]|| 2 3 as an example of the way we prove the above estimates. We have 



||PDAqz]M)|| 2 < \\h(z)\\ 2 s ( sup|fc|V^e-^) [ds 

\ |fe|<l,T>0 / JO 



t - s 



|D/;(z)||o _ { sup e ^) fds (1 + S)_ 

|fe|<l,r>0 / Ji 



^^wii fy^+TTTT^ r-^)<ciizn 2 (i+trt ,2.15, 

for all t > 0, which shows that ||PD7V[z]|| 2 3 < C||z|| 2 . All other estimates in ( 12.14b can be done similarly; we postpone 
their proof to Section|6]below. 

Finally, we note that the Lipschitz-type estimate in ( 12.131 > can be obtained in the same manner, mutatis mutandis, due to the 
similarity between (|2.91 l and ( 12.101 ) with ( 12.81 ): we omit the details. ■ 

We can now turn to the question of the asymptotic structure of the solutions of (11.11 ) provided by Theorem [6] Note that 
already if we wanted to prove that e^*Zo satisfies 'Gaussian asymptotics' we would need more localization properties on 
zo than those provided by the So-topology. It will turn out to be sufficient to require zo G Bq n L 2 (IR, x m dx) for (some) 
to > 2. We now prove that this requirement is forward invariant under the flow of dl.lt : 
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Proposition? Let p m (x) = \x\ m and define 

V rn = | z e B such that |z | + |||O m Zo||2 < oo | . 

Ifzg £ T> m and \z$\ < £o such that Theorem\6\holds, then the corresponding solution z(t) of ( 17.71 ) satisfies z(t) £ T> m for 
all finite t > 0. Furthermore, there holds \z(t)\ < (1 + <$)eo /or some (small) constant S. 



Proof. Note first that by Theorem|6j |z(i)| < ||z|| < (1 + <J)eo since zo 6 Bo and |zo| < to. Then, fix m £ N, m > 1. The 
proof of Theorem [6] can easily be adapted to show that dl.lt is locally (in time) well posed in T> m . Global existence then 
follows from the fact that the quantity 



1 1 f°° 

N(t) = -\\p m z(;t)\\ 2 =2j dx \x\ m (a(x,t) a + b(x,t) 2 ) 



grows at most exponentially as t — > oo. Namely, we have 

d t N(t) = fdx \x\ m (d x (ab) + 2b8 2 x b + bdjf(a, b)d x b + g(a, b))) 

coo 

dx m|a;r" 1 sign(a;)^(a + g(a, b)) + (2 + /(a, b))bd x b 
dx\x\ m (d x b) 2 (2 + f(a,b)) 

-oo J 

< fdx ((m - I)" 1 ' 1 + \x\ m ) b(a + g(a,b)) + (2 + f(a,b))bd x b 

J — OO 

- [dx\x\ m (d x b) 2 (2 + f(a,b)) 

J — oo 

< [dx ((m - I)™" 1 + |x| w Y|6(a + 9{a, b))\ + 2^|2 + /(a, 6)|6 2 ) 
<C 1 (m,e ) + C 2 (e )N(t) , 



due to the estimates ||/(a, &)||oo < C £ o "C 2 and 



<C*e . 



3 Asymptotic structure - Proof of Theorem 21 

We can now state our main result on the asymptotic structure of solutions of dl.ll ) in a definitive manner: 

Theorem 8 Let T> m be as in Proposition\7\with m > 2, let Zq £ D m with |zo | < eo such that Theorem\6\holds and define 

u(x, i) — a(x — t,t) + b(x — t, t) and v(x, t) = a(x +t,t) — b(x + t, t) 

for the corresponding solution z(t) — (a(t),b(t)) of ( I7.7D . Then there exist functions {gn}n=o (independent of zq) and 
constants Cm, {dn)n=i determined by z$ such that 

W ~ (3-1) 
wfof) = - 7 =L= ff -(-^_) + £ * ) + jR f (a;,*) , 
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where the remainders and R~] satisfy the estimates 

sup(l + t)i-^||i?f Ui „ } (-,t)|| £2(R) <C N 

"° , t (3.2) 
Bup(l + t)* ^\\d x R( u . v} (-,t)\\ L 2 < C N . 
t>o v ' 

Furthermore, for n > 1, the functions satisfy <?„ (z) ~ 

\z\- 1+2 ~ n asz — > ±00. 

Remark 9 As will be apparent from the proof of Theorem^ any hyperbolic-parabolic system of the form 

d t z + f(z) x = (B(z)z x ) x 

with admissible nonlinearities in the sense of (the natural extension of) Definition\l\gives rise to solutions having the same 
asymptotic structure as those of the p-system as long as the following two conditions are satisfied: 

1. There exist two matrices S and A with S non-singular and A diagonal having eigenvalues of multiplicity I for which 
Se^ f w e^ ot S in the sense of LemmaU9\(see Section^, where Lq — d 2 x + Ad x and L — B(0)d% — f'(0)d x . 

2. The Cauchy problem with initial condition in the corresponding functional space (the natural extension of Bo to the 
problem considered) is well posed and satisfies the analogues of Theorem\6\and Proposition^ 



We now briefly comment on the above assumptions for specific systems such as the 'full gas dynamics' and the MHD 
system. The intertwining property of itemQ]above is proved in |4| for quite general systems, though not in exactly the same 
topology as that used in Lemma [T9l As for item|2l local well-posedness for initial data in Bq is certainly not an issue, the 
only difficulty is to prove that the various norms of Definition[5]exhibit 'parabolic-like' decay as t — > 00. This is very likely 
to hold, particularly for systems satisfying itemQ] 

While the variables (a, b) are adapted to the study of the Cauchy problem because of the inherent asymmetry of spatial 
regularity in dl.ll ), they are not the best framework for studying the asymptotic structure of the solutions to dl.lt . It turns 
out to be more convenient to change variables to quantities that move along the characteristics. We thus define 



\Sz(x, t) 



fu(x,t)\ = (T- 1 0Yl 1 \Mx,t)\ = (T~ l OY 
\v(x,t)J-\0 T\l -lJ\b{x,t)J ~\ 7/ 

where T is the translation operator defined by 

(Tf){x,t) = f(x + t,t) or equivalent^ by ff(k,t) = e lkt f(k,t) . (3.3) 

Note in passing that 

a(x, t) = —(^u{x + t,t) + v(x — t, t)\ and b(x, t) — — (u[x + t,t) — v(x — t, t) \ . 
We then use the fact that z satisfies the integral equation 

5z W= 5cLtzo + J (^5e^)(^ /i( z(s)) ) 

Sz + Jd S e^s( d J {z{s)) )+n[z}(t), (3.4) 



e L °'' 



where 
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To justify the notation, which suggests that 1Z is a remainder term, we will prove in Section [7] that lZ[z] = (1Z U [z] , 1Z V [z]) 
satisfies the improved decay rates 

l|ft{u,i>}Nlkr + ll D ^MNII 2 ,r ^ Ce ° » ( 3 - 5 ) 

because of the intertwining relation <Se^* w e^°*5 (see Lemma |T9T> and the fact that h(z) — (?o( z ) + h.o.t.. 

Recalling that g is quadratic (cf DefinitionQ]), we will write 

g (z) = c+(a + bf - c_(a - 6) 2 + c 3 (a + 6) (a - b) 
= c+{Tuf - c-{T- l v) 2 + c^Tu^T^v) 

for z = (a, 6). We thus find from ( 13.41 ) that u and v satisfy 

u{t) = e a «*(a + 6o) + jjds e d " {t ' s) (c + u{sf - c_r- 2 w(s) 2 ) 

+ r- 1 ^ M [z](t) + C3 a :c ^d S e^( t - s )r- 1 ((r M ( s ))(r- 1 u ( s ))) , (3.6) 

v{t) = e 9 '*(a - 6 ) + 9, Jds e a -('- fl >(c_i;(s) 2 - c + T 2 u(s) 2 ) 

+ r^[z](i) - C3 d x J As e 9 '^ s )r((T U ( s ))(r- 1 w ( s ))) . (3.7) 

Note that, but for the presence of the second lines in (13.61 ) and ( 13.71 ). these expressions are precisely Duhamel's formula for 
the solution of the model problem jl.71 ). written in terms of u = T~ lr u + and v — Tu^. The next step is to write 



u = + R u = uq + u\ + R u and v = v* + R v = Vq + v % + , 
considering R^ and R„ as new 'unknowns' and 

1 N 1 

«o(s,t) = -7=^(7^) , ui(x,*) = J] 7 u - 1 

1 N 1 

«o(a;,i) = -==5o and «i(*»*) = 2Z- 7i — d n9n(-wr= t ) 

VI + t Vi+t ^(l + t) 1- " - " 



(3.8) 



2^+1 



for some coefficients {d^}^ =1 and functions }^L to be determined later. 
We now use 

W 2 = (tt — + U*) + U 2 = R U {U + U*) + U 2 + 2lioMl + Uq , 

« 2 = (v - v*)(v + u*) + = i?f (« + u*) + w 2 + 2<u ui + w 2 , 

(T M ) ( r-,) = (Ti?f)r- 1 (^)+ (r I «f ) r(^) + cnaT-V) 



Since 



</+(:r) = u (x,0) , Ui(x,0) = ^d+g+{x) 

71=1 
N 

5q(x) = v o (x,0) and i> x (a;, 0) = ^ d~g~(x) 
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we find that J2„ and satisfy 

R»(t)=e a * t (a + b Q ~g+) 



s^uoiO) + c+d x / ds e a ' ( *- s) M (s) 



u (t) 



■ d * t u 1 (0) + 2c+d x /dse a - (i - s) u (s)iii(s) - u x (t) 
Jo 

r r t -I JV 

c_ a a ds e d ^T- 2 ((v ( S ) 2 + 2vo( S )v 1 (s))) - £ e^*d+ S + 
L Jo v y J „ =1 



(3.9) 



i^(t)=e^ t (o -6o-^) 



e^voiO) + c^d x I dse d * (t - s) v (s) 2 



vo(t) 



^*wi(0) + 2c-d x / ds e 9 *(*-%o(s)«i(s) 



d x /d S e 9 ^ i - s )T 2 (( Uo (s) 2 + 2uo(s> 1 ( S ))) -]Te a ^ 5 

«/ J ~ i 



^[z,R iV ](i) + r^[z](t) , 



where 



■R u [z,R N }{t) = c + E [^i, u + h 3tU ](t) - c_E_ 2 [/n,„ + fe^K*) + c 3 E_i[/i 2 + MO) 
ft„[z, R N ](t) = c_E [/i M + ft3 il( ](t) - c + E 2 [/n, u + ft 3 ,«](i) - c 3 E![/i 2 + 7i 4 ](*) , 
withR Ar = (i?^,i?f), 



Ea[h}(t) = d x ds e d ' {t ~ s} T a h{s) and 
Jo 



h ltU = R%{u + u*), h 3 , u = ui, h 2 = {TR£)T 



N\n--i( v + v* 



(T~ 1 R®)T\ 



u + u* 



(3.10) 



hi, v = R v (v + v*) , h,3,v = v x , hi = (Ttt*)(T u*) . 

Note that we can write (fJU) and d37TOb as H N = R N }. If we now consider z fixed, we can interpret R^ = J-[z, R N ] 
as an equation for R w which can be solved via a contraction mapping argument. Namely, we will prove that if ||z|| < Ceo, 
H N i ► !F[z, R N ] defines a contraction map inside the ball 



\R N 



I "J _ 4 • r llJ.-f-' II"'' l-'-T ' " Ikl-e 



pJV| 
it,, 



< c 



(3.11) 



fore = 2 N 2 , provided {g^}^ =0 and {d^}^ =1 are appropriately chosen. 



Basically, we will choose uo, fo, ui and «i in such a way that the second and third lines of ( 13.91 ) and ( 13.10b vanish. Note 
that if, for instance, we set the second, respectively third lines of ( 13.91 ) and (13.101 ) equal to zero, the resulting equalities are 
nothing but Duhamel's formulae for Burger's equations for u and v , respectively for linearized Burger's equations for u x 
and v\. Properties of solutions to these types of equations are studied in detail in Section|4]below. 

Once Mo, vq, u% and vi are fixed, the time convolutions in the fourth lines of ( 13.9b and (13.10b can then be viewed as the 
solution of inhomogeneous heat equations with very specific inhomogeneous terms. Properties of solutions to this type of 
equations are studied in detail in Section[5]below. 



Assuming all results of Section|4]and|5] we now explain how to proceed to prove that T\z, R N ] defines a contraction map. 
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Obviously, the requirement on {gn}n=o an d {dn}n=i i s tnat tne m " st f° ur li nes m ( 13.91 ) and ( 13.10b satisfy ( 13.111 ). This is 
achieved in the following way: 

1 . The first line of (13.91 ), respectively of (13.101 ) satisfies ( 13.1 II ) for any such that the total mass of g^ is equal to that 
of ao ± bo, provided ao ± bo and g$ satisfy || x 2 (ao ± £>o ) 1 1 2 < 00 an d II a;2 ffo Ih < °°- This fixes the total mass of 
g . Note also that we need the estimate || x 2 (ao ± £»o ) 1 1 2 < 00. There is no smallness assumption here, which is to 
be expected since generically || x 2 (a(-, t) ±b(-,t))\\2 will grow as t — > 00. Note on the other hand that Proposition|7] 
shows that || x 2 (a(-, t) ± b(-, t))\\2 remains finite for all t < 00, so requiring || x 2 (ao ± fc>o ) 1 1 2 < 00 is acceptable. 

2. We can set the second lines in ( 13.9b and ( 13.10b equal to zero by picking for uq and vo any solution of Burger's 
equations 

d t u = dliiQ + c + d x (u ) 2 and d t v = d%v + c_d x (v ) 2 

(or of the corresponding heat equations if either c+ or c_ happen to be zero). In Proposition [12] we will prove that 
there exist unique functions uq and vq of the form given in ( 13.8b that satisfy the conditions of item[TJabove (total mass 
and decay properties). This uniquely determines uq and vq. 

3. We can also set the third lines in ( 13. 9b and (13.10b equal to zero, by picking any solutions u\ and v\ of linearized 
Burger's equations 

d t U\ = d 2 ui + 2c + d x (u ui) and d t vi — d 2 vi + 2c_d x (v Vi) . (3.12) 

In PropositionfTSl we will also prove that there is a choice of functions sucn tnat Ul an< ^ Vl m <EU satisfy 

(13. 12b for any choice of the coefficients {d^}^ =1 . Furthermore, in Proposition IT2l we will prove that the choice of 
functions can be made in such a way that g^ (x) have Gaussian tails as x — > q=oo and algebraic tails as x — > ±00. 
This actually completely determines g^ (x) up to multiplicative constants (this last indeterminacy will be removed 
when the coefficients {d^ }^ =1 are fixed). 

4. We then further decompose the terms involving in the fourth lines in (13. 9b and (13. 10b as a^r (x) — f n (^fx)+R^(x). 
The definition and properties of f n (x) are given in Lemma [TOl In particular, in Proposition [12] we will prove that 
Rn{x) have zero total mass and Gaussian tails as |a;| — > 00, which implies that e^R^ also satisfy ( 13. lib . 

5. Finally, in Section [5] we will prove that the time convolution part of the fourth lines in (13. 9b and ( |3.10b can be split 
into linear combinations of e 9 ** fn(Tx) with n = 1 . . . N + 1 plus a remainder that satisfies ( 13.1 lb . The coefficients 
{ d-n } n= 1 can tnen be set recursively by requiring that all the terms with n — 1 ... N coming from the time convolution 
are canceled by those coming from item|4]above. This can always be done because the coefficient of e^'/ m (=Fx) in 
the time convolution part of the fourth lines in ( 13. 9b and ( 13.1 0b depends only on g^ if m = 1 and on d^ l _ 1 if m > 1. 
The only term that cannot be set to zero is the last term in the linear combination (the one with n = N + 1), which is 
the one that 'drives' the equations and fixes e = 2~ N ~ 2 . 

The procedure outlined in[TJ[5]takes care of the first four lines in (13. 9b and (13.10b . We will then prove in Section [7] that the 
terms 7?-{ Ui „} [z] satisfy ( 13.11b and that 

1 1 
J2 \\V a n { u, v }[z, R N ]\kl +§ -e < Ce £ ||D«R"|| a>m _ e + C , (3.13) 

a=0 a=0 
1 1 

£||D a (^ { „, l)} [z,Rf]-^ K , )} [ Z ,Rf])|| 2i 3 +f _ e <C'eo£||D-(Rf - R*)|| 2if+f _ e . (3.14) 

a=a a=0 

This finally proves that J-[z, JH N ] defines a contraction map and that the solution of R N = J-[z, H N ] satisfies ( 13. lib , which 
completes the proof of Theorems [2] and [8] 
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4 Burger's type equations 



In this section, we consider particular solutions of Burger's type equations 

d t u =dlu + 7<9 x Mq 
d t u„ = dlu^ + 2jd x (u Un) 



of the form 



u o(x,t) = -2—g (-^—) and u±(x,t) 



(l+t) ^ 



< _g±( *) 



(4.1) 
(4.2) 



(4.3) 



/oo /> OO 

dx w (x, t) = dx go(x) small enough, there is a unique choice of g and (^j 
-oo J — oo 

such that (x) — f n {T x ) + Rn i x )> where 



fn(z) = / d£ 



(4.4) 



and has zero mean and Gaussian tails as |x| — > oo. In particular, g£ (x) decays algebraically as x — > ±oo, as is apparent 
from (l4~4t , 

Before proceeding to our study of ( 14. lj and (14.21 1, we prove key properties of the functions /„. 



Lemma 10 Fix 1 < n < oo. The function f n is the unique solution of 

d 2 J n (z) + \zd z f n {z) + (1 - ^r)fn(z) = , With 

UO) = 2^r(i±|l!l) ant/ lim z" 1+ *cT/„(z) < oo . 

/oo 
dz = and there exists a constant C(n) such that 

- oo 

2 

sup £ P£_ m , 1+m _^(«)|flT(*/»(*) + WzUz))\ < C(n) 

sup ^ P^-i- ro ,2+m-^-(*)I^T/n(*)l < 



Pp,<f(z) = 



(l + z 2 )§eT i/;z>0 
(l + z 2 )i ;/z < 



(4.5) 



(4.6) 



Proof. We first note that /„ can be written as 



fn(z) = / d£ 



2 / d^^^dtfe-- 



(g + z)e ~ 

io i 1 -^ Jo 

This shows that /„ solves i4.5\ since, defining Cf = d^f + \zd z j + (1 - ^rr)f' we nn d 



Cfn(z) 



[e^a|(e- l£ ^)- 5 ^ T (-2)^- 1 ^(e-- 



(4.7) 



13 



Obviously, f n (z) is finite for all finite z, so we only need to prove that /„ satisfies the correct decay properties as \z\ — ► oo 
so that (14.61 ) holds. It is apparent from ( 14.41 ) that f n decays like a (modified) Gaussian as z — ► oo and algebraically as 

z 2 

z — > — oo. Furthermore, substituting /(z) = C|z| Pl and /(z) = C|z| P2 e~T into Cf = shows that the only decay rates 
compatible with Cf = are p\ = — 2 + and p 2 = 1 — ^r- 

We now complete the proof of the decay estimates (14.61 ). Let F n . m (£,z) = <9™((£ + 2:)c _<?+ 4 ) ) and G n , m (£,z) = 
8?{zF n (S,z) + 2d,F n (Z,z)). 

We first consider the case z > and note that F n , m and G n , m satisfy 

\F n , m {£,z)\ < \F n<m (0,z)\ and |G n ,„,(£, < \G n , m (0, z)\ 
for all £ > if z > zq for some zq large enough. We thus get, e.g. 



!/»(*)! = 



<\F n;O (0,z)\ / dC^-i + z 1 -^ / d£|i^, (£, 2 )|<C^-^e-^ 



The estimates on |<9™(z/„(z) + 29 Z / T j(z))| and \dl +m f n (z)\ when z > and m > 1 can be done in exactly the same way; 
hence we omit the details. 

We now consider the case z < and note that F n , m and G n , m satisfy 

\F n>m (Z,z)\ < \F nim (-§,z)\ and |G n , m (£, z)\ < |G„, m (-f , z)\ 
for all < £ < — | if z < — z for some z large enough. We thus find (integrating by parts in the second integral below) 

|/„(*)|= /dCF^o^z)^" 1 



<|F r , (-f,z)| / d^" 1 



< C|«|* _1 e - fe + 2(1 - 4 ) / d£ e < c^l^ 



Since the remaining estimates can again be done in exactly the same way, we omit the details. It only remains to show that 
f n {z) has zero total mass. This follows from 

SO /-OO 

dz f n {z) = (§- / dz Cf n (z) = , 



since 9f / n , zd z f n and /„ are all integrable over R. ■ 

Remark 11 Using the representation \4.7\ , splitting the integration interval into [0, 2~~z) and [2 _ ~2,oo), integrating by 
parts and letting n — ► 00, one can prove that 

lim 2-"/„(z) = ze-T ■ 

n — >oo 

which shows that the constant C (n) in \4.6\ grows at most like 2 n . 



We can now study in detail the solutions of ( 14.11 ) and ( 14.21 ) that are of the form ( 14.3 



Proposition 12 Fix 1 < n < 00. For all a, 7 S R with \aj\ small enough, there exist unique functions uq and of the 
form ( 14.51 ) that solve ( 14.71 ) anc/ ( I4.2D . w!//z go satisfying 



dz go(z) = a 



x - e 4 



,2W 



|^*Sb(z)| < C|a| 
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and with g^(z) — f 7 i(T z ) + Rn ( z )> where satisfy 

/oo 3 2^- 

dzR±(z) = and sup ]T — ==^— — l^nWI < C|a 7 | . 

POO 

Proof. The (unique) solution of d4.ll > of the form Uq(x, t) = , g ( * ) satisfying / dz go(z) = a is given by 

_ tanh(^)e~^ 

~ 7 A(l+tanh(^)erf(f)) ' 

In particular, we have 

tj-^=^\dT9o( Z )\<C\a\. (4.8) 

We next note that substituting ( 14.3b into ( 14.21 ) gives 

= d 2 5 ±(z) + \zd,cfc{z) + (1 - 3^)<£(*) + 2ld z (g a (z)g±(z)) 
= £g±(z)+2jd z (uo(z)g±(z)). (4.9) 
We formally have (using integration by parts) 

dz gt{z) = (| - / dz £ ff ±(z) + 2 7 d z Mz) 5 ±(z)) = , (4.10) 

-oo J — - oo 

which shows that g,^ have zero total mass, provided the formal manipulations above are justified, i.e. provided g^ and its 
derivatives decay fast enough so that the integrals are convergent. 

As is easily seen, f n {z) and f n (—z) are two linearly independent solutions of Cf = 0, whose general solution can thus 
be written as c\f n (z) + cij n {—z). Using the variation of constants formula, we get that the solution of d4.9b satisfies the 
integral equation 

ate) = uz)(ci + 2,jk i^o^Mmy fni _ z) (4 27 j^ wyi ) , 

where the Wronskian V^(z) is given by W(z) = f n {z)d z f n (—z) — f n (—z)d z f n (z) and c x and are free parameters. 

Note that W(z) satisfies d z W(z) — —^W(z) and hence W(z) = W(0)e~~ for some W(0) ^ 0. We now set cf and cf 
in such a way that (after integration by parts), we have 

g±(z) = f n (Tz) + R[g±](z), (4.11) 

R[9i](z) = ffoy/nW f d£ e^(e/„(-0 + 2^/ n (-0)5o«)5^(0 

^ oo 

/OO 2 
^eV(e/ n (e)+2^/ n (o)sb(o^(0- 

Using Lemma[lO]and d4.8l ), it is then easy to show that for |cry| small enough, d4.1 It defines a contraction map in the norm 

i/i 2 _^su P (yrT^) 2 -^i/(z)|. 

Namely, we have the improved decay rates 

l si 

sup , — 4u +m _ x iCTffnK*)! < C|«7l |^| 2 _^ • 
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This shows that d4.1 U has a (locally) unique solution among functions with |/| 2 _^ < cq if |cry| is small enough. In 
particular, there holds 

l si 

zeR m^0 W 1 + Z ) ^ 

from which we deduce, using again ( 14.111 ) and LemmafTOl that |Dg r f |3_^. < c\ and thus 

sup fl , \dlR[gt]{z)\ < C\crf\ . 

zgr (vi + z 2 y~^ 



Iterating this procedure shows that |D m g„ | 2+m _^_ < c m and that 



3 Ml 

e 4 



sup ^ , \8?R\£]{z)\ < C\a 7 \ 

as claimed. In turn, this proves that the formal manipulations in d4.101 i are justified, so that the functions (z) have zero 
total mass, which shows that the remainders R[g^](z) have zero total mass as claimed since R[g^](z) — g^(z) — /„(±z) 
and since both g^ (z) and f n (z) have zero total mass. ■ 



5 Inhomogeneous heat equations 

In this section, we consider solutions of inhomogeneous heat equations of the form 

ftu = ^u + e ((l + i)*-*/(^2)) , u(a:,0) = 0. 
where / is a regular function having Gaussian decay at infinity. Solutions of ( 15. U satisfy 

Theorem 13 Let 1 < n < oo, a = ±1, E(x) = c~, M(/) = dz f(z) and 

J — oo 

i 1 roo i 



«n(M) = ° - i 2 /n(vB) ^ /»W - / d£ / u ; ■ (5.2) 

(l+t) ^TT V4tt _/ z (^-z) 1 -^ 



(5.1) 



The solution u of ( 15. 1 i satisfies 

2 

||u-M(/)«„|| 2iJ * + ||D(«-M(/)u n )|| 2)r < C V HSD-/IU , (5.3) 

m=0 

/or / 5mc/j f/iflf f/;e r./z.i. of ( 15.31 ) is finite. 

Remark 14 Atofe f/iaf while u — » M(/)u n as i — > oo /n f/ie Sobolev norm ( 15. 3D , i/ ofoes nof afo so /n spatially weighted 
norms such as i 2 (R, x 2 da;), as u„ /zas infinite spatial moments for all times, while all moments of u are bounded for finite 
time. 

Proof. We first define 

F(0 = f dz (f(z) - M(f) ^1) with M(/) = /°dz f(z) (5.4) 

J — oo J J — OO 
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and note that F satisfies 

2 2 2 

l|D 3 ^l|i + E \\p^ m F\\i + ]T \\V m F\\ 2 <CJ2 l|2D ro /IU , (5.5) 

m— m—1 rn—0 

where = \/l + x 2 . Namely, we first note that \\pF\\i < \\F\\ 2 + \\F"\\ 2 and F(k) = (ik)- 1 (f(k)-f(0)e- k2 ). Then, 
since IjS/Hoo < oo implies that / is analytic, F is regular near k = 0. The proof of ((53) now follows from elementary 
arguments. 

We finally note that it follows from (15. 4t that 

(1 + t)^-if(^^Si^= M(/) 2 e - ( Vi+t) 2 + (l + t)^-!^^^^) . (5.6) 



=A(s,t) =da,B(x,t) 

The proof of d5.3t is then completed by considering separately the solutions of heat equations with inhomogeneous terms 
given by d x A(x, t) and d 2 B(x, t). This is done in Propositions [T5land[T6lbelow. ■ 



Proposition 15 Let <r = ±1, 1 < n < oo, and let u n be defined as in ( 15. 2D . The solution u of 

d t u = d 2 x u + d x A , u(x, 0) = , (5.7) 

with A defined in ( 15.61 ) satisfies 

\\u - u n \\ 2> 3 + ||D(« - Wn)|| 2 ,| < C . (5.8) 

Proof. The solution of d5.7| ) is given by 



r e 4 <*- s > e 4 < 1+s > 

u(a;,i) = d x ds dy—=== — — -3— - . (5.9) 

Jo J-oo V 47r (* ~ s ) v4tt(1 + s)2-2^ 

To motivate our result, we note that performing the ^-integration and changing variables from s to £ = 2 ^~^ £ in ( 15.9b leads 
to 

1 ~ I ~ f 2 

lim(l + i) 1 -^Vr U (- ( rzVm,<) = lim 22^± / ^ d£ * c ~ 4 ; i = ct2 ""^ /»(z) . 

More formally, taking the Fourier transform of ( 15.91 ) gives 



\{k,t)=ikc- k ( 1+t > ds 

Jo (l + s) 1 



We now use that 



t Jlikas 

ds ; I ds 



t Jlikas 



< C{n) , 



ds— — = \k\-^(e(ak)J n (\k\t) + e(-ak)J n (\k\t) 
Jo s 2 " v 

where 9(k) is the Heaviside step function and we defined 



J n {z) = ds — — - 
Jo s 271 
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for z > 0. This function satisfies 

supz 1_ ^| J n (z) - J n ,oo\ < 7: for Jn.oo = lim J n {z) . 
z>0 £ z ^°° 

Now define 

u^(k,t) = ike- k2 ^\k\-^(0(ak)J 7hOD +e(-ak)X^^j . (5.10) 

We have 

\u(k,t)-u^(k,t)\ < {C(n)\k\ + i- 1 +^) e - fc2 ( 1+t ) < {C(n)\k\+r^)e- k2{1+t) , (5.11) 

from which ( 15 ,8b follows by direct integration. We complete the proof by showing that the inverse Fourier transform of the 
function u^(k, t) defined in d5.10t satisfies 

1 />oo p — S_ 

Un (x,t) = ° Uf=) for /«(*) = / d£ -i— ; . (5.12) 

This follows easily from the fact that 

t£(M) = (l+i)^ + ^ r %(tv / lTt,0) , 

and that, since 



1-- 



1 5 



we get 

£^L /n( _ afc) = 2 -^zA:e- fe / d£ * = ife e - fe / d£ - , 

= ike- k2 \k\-^[9(ka)J ni0O + 6{-ka) J^j= u^(k, 0) 

as claimed. ■ 

Proposition 16 Lef <r = ±l, l<n<oo and p(x) — VT+ x 2 . 77ze solution u of 

d t u = dlu + d 2 x B , u(x,0)=0, (5.13) 

w/f/i _B defined in ( 15. 61 ) satisfies 



Hkf* + ||Du|| 2)r < C M|D 3 f ||! + J] ||pD ro f ||! + J] ||D m F|| 2 j (5.14) 

^ m— m— 1 ' 

/or a// F for which the r.h.s. of \5.14\ is finite. 

Proof. We first note that the Fourier transform of u is given by 

/•* 

u{k,t) = -k 2 ds c-^its^kas P(ky/T+!)(1 + a)*"* , 

Jo 

which implies 

'* ds 



||(1-Q)«|| 2 3 +||(1-Q)Dm 



| 2 ,| <C(||DF|| 2 + ||D 2 ^|| 2 ) sup / 

4 v '0<i<lVo 
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Here Q is again defined as the characteristic function for t > 1. Next, integrating by parts, we find 

u(k, t) = - - ^(kVTTt)e^ + ^ 

where N(k, t) = £ [ds ^it- s) -^s( e + d s ) ( . 

2ct Jo V /\(H-s)3-2^y 



We then note that 

||u - iV|| 2 , s + ||D(« - A0|| 2)J < c(\\F\\i + \\DFh + l|D 2 F|| 2 ) , 
and that, defining G(k) = \d k F{k), we have N[k, t) = N {k, t) + N^k, t) + N 2 (k, t), where 

iVo(M) = - fdse-^-^ (i^HH. 
2^ J q V(l + s) 



1 _ 1 / 3 

2 2"" 



2ct Jo Vfl + s) 1 "^/ 



N 2 {k,t) 



2ai 2 - 5 ;y dse 1(1+^-*/ 



The procedure is now similar to that outlined in the proof of Theorem [6j split the integration intervals into [0, |] and [|, t] 
and distribute the derivatives (fc-factors) either on the functions F and G, or on the Gaussian. Introducing the notation 



^n(t) s /d s li+i^!i+ / ds ii±£)_, (5.i5) 



-Q2 



we then find that 



|QD«iVo|| 2 , f+f < C(\\F\\i + ||D 2 +«F|| 1 ) S upti+f S x [I+Jf ] (t) 



t>i 



°iv 1 || 2ii+9 <a(||G|| 1 + ||D 1 + 'G|| 1 ) 8 u P t*+* BitltyiW, 

4 2 t>l L 4' 1 +2 

^i kl+r < cmi + iid^iio-p^I^ b x [}2J] (*) 

for a = 0, 1. The proof is completed by a straightforward application of Lemma [18] below, where we consider generaliza- 
tions of the function B\ in ( I5.151 l (see Definition[T7lbelow). since those will occur later on in Sections|6]and|7] ■ 



6 Proof of Theorem [6], continued 

In view of the estimates ( 12.61 > and ( 12. 81 > on e^' and /i, respectively, the estimates needed to conclude the proof of Theorem 
|6]will naturally involve the functions Bq and B which are defined as follows: 



Definition 17 We define 

B [q)(t) = [ds 
Jo 



o V^^(l + s) 9 

](t)=/j 7 £ + ^ , + /ds (1 + a) ~ 92 ln(2 + S ( 3 . (6.1) 
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These functions satisfy the following estimates: 



Lemma 18 Let < p 2 < 1, < r 2 < 1 — p 2 , Pi) Sij <fe, ?"i > ana? r^, £ {0, 1}. There exists a constant C such that for 
allt>0 there holds 



B [«i](t)<C(l+t)-« 



i 



i/ < pi < 1 



w/iere /3 = min(p! + min^ — 1, 0) + r%,p2 + q 2 + r 2 — 1). — m a x (^qi,ij <5p 2 +r 2 .i + r 3) an d is the Kronecker delta. 
Furthermore, since 

the estimate in ( 16. 2D applies for B\ as well. 

Proof. The proof follows immediately from 

^ r w x t ds 1 /"f e~* 

Bobi]( ^ c ^y o ^ + (iTTFy ds 7i' 

Lp2,92,r 3 ,r 3 J W - (|)Pi(* y (l + a)9i (* + 1)92 _/ S P2 (1 + s )r 2 

and straightforward integrations. ■ 

We can now complete the proof of Theorem[6] 

Proof of Theorem® continued. 

First, we recall that our goal is to prove that the map M defined by 

satisfies ||A/"[z]|| < C for all z e S with ||z|| = 1. We have already proved that ||PD7V[z] || 2 a < C. The other necessary 
estimates are done as follows: 



i i 



\\^]\\oo,o< C sup B 1 [fl](t)<C, 

t>0 2 ' 2 

\WM\\ 2i <c , su P (i + i)3s 1 [! , !](t) <c, 

||PDA/"[z]|| 2 3 < Csup(l + t)3fli[i! l(t) < C, 

' 4 t>o 2,4 

||(l-P)DAA[z]|| 2 3 <sup(l+*)*B [f](t)<C, 
4 t>o 

||(l-Q)PD 2 AA[z] 2 || 2 s * <C||(l-Q)PD^[z] 2 || 2 3 <C||PDAA[z] 2 || 2 3 <C, (6.4) 

'4 '4 '4 

||QPD 2 AA[z] 2 || 2jr <Csu P £±a|i?[ff ( ; ](t)<C, (6.5) 



||(l-P)D 2 AA[z] 2 || 2r <sup(l + t) 4 B [|](t) <C. (6.6) 
4 t>o 
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In (El*, we used the obvious estimates ||PD/|| 2 < ||P/|| 2 and || (1 -Q)f\\ 2 , P < 2 p - q \\ (1 -Q)/|| 2 , g if q < P, while in O), 
we made use of sup \k\\/l + te~ k * < 1, and finally in ( 16.6b we used sup fceR |fc|(l + fc 2 ) - ^ = 1, Incidentally, (16.61 ) is 

|fe|<l,t>0 

the only place in the above estimates where the (crucial) presence of the extra factor (1 + fc 2 )~ 5 in the second component 
of the r.h.s. of ( 12.61 ) is used. This concludes the proof of Theorem[6] ■ 



7 Remainder estimates 



We now make precise the sense in which the semigroup e^* is close to that of (12. 3K whose Fourier transform is given by 

-k 2 t+ikt q 



Xot _/e 







-k 2 t-ikt 



(7.1) 



Lemma 19 Lef P fee f/ie Fourier multiplier with the characteristic function on [—1, 1], one/ Zef e^* resp. e^°* fee as /n ( 12.21 ) 
resp. ( 17. 7 1 ) ant/ 5 fee as z'# ( I2.4D . 77ien one has the estimates 



sup vi + te 2 
t>o,fceR 



P5e i4 - e Lot S 



where (PSe Lt - e Lot S)i, 3 denotes the {i,j)-entry in the matrix PSe" - e u,t S. 



1,3 



<C, 



(7.2) 



Proof. The proof follows by considering separately \k\ < 1 and \ k\ > 1. We first rewrite 



P5e Lt - e Lot S = Pf 5e Lt - e L °*<S)+ (1 - P)e Lo *5 . 



We then have 



sup \f \ + te 2 
oo.fceR 



(1 - F)e Lot S 



< sup y/l+te V < C 
t>o,|fc|>i 



For |fe| < 1, we first compute 



° 1 I e -ifct _ e -ikt | > 



L , ; , t cos(fctA) + i sin(HA) cos(fciA) + ±^ z sin(HA) 



<Se Lt = e" 



^os^tA) - i^-i sin(feiA) -(cos(fciA) - sin(fetA)) 



where we recall that A = yl — A; 2 . We next note that 

P| sin(fctA) - sin(fci)| + P| cos(fciA) - cos(fct)| < P| cos(fci(A - 1)) - 1| +P| sin(fc<(A - 1))| 

< Plv 7 ! - fc 2 - 1| \k\t < ¥>\k\ 3 t , 



i - l)sin(fctA)| < Plv 7 ! - k 2 - 1| < P|fc| 3 t . 



The proof is completed noting that 



sup t^\k\ n e~— <C(n) 
l*l<M>0 
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for any (finite) < m < n. ■ 

We are now in a position to prove that the remainder 







)zo + 


[ds 




Jo 



\ _ e L (t-s) s ( 
h(z(s))J \d x g (z(s)) 



satisfies improved estimates as stated in (13.51 >: 



Theorem 20 Let eo be again the ( small) constant provided by Theorem® Then for all Zq G Bq with \zq \ < eo, the solution 
z of di.il ) satisfies 

||ftMlla,!* + l|Drc[*]||2,f* <Ce . (7.3) 

Proof. We first note that 

(Se Lt - e Lot s)z =(sPe Lt ~ c Lot <s)z + 5(1 - P)c Lt z = Li[z ](*) + L 2 [z ](t) , 

and then use the fact that by Lemma [19] we have 

||D Q Li[z ]|| 2 < Csup(l + t)Ht min(||D a zo|| 2 , H-^||fo|U)< C|z | 
' 4 2 t>o v ' 

for a = 0, 1 and finally 

||L 2 [z ]|| 2 a + ||DL 2 [z ]|| 2 6 < C(||z || 2 + ||Dzo|| 2 )sup(l + t)U~t < C\z \ . 

t>o 

This proves 

||(<Sc Lt - e L o*s)z || 2!f + ||D(5e L * - e Lot s)z || 24 < C|z | 
for all zo £ Bo. We then show that 

\\K{z](t) ~(se U - e L ° t s)z \\ 2 , r + \\B(n[z](t) -(Se Lt - e L ^s)z )\\ 2 ^ < C\\z\\ 2 
for all z G B. We only need to prove the estimates for ||z|| = 1. We first decompose 

K[z](t) -(Se Lt - e Lot s)z - SA/i [z] (t) + SM 2 [z] (t) + -A/" 3 [z](i) , (7.4) 



where 





d x h(z(s))J ' 



M[z](t) = (l-F) /dse L( *- 5 
Jo 

W 2 [z](t)=P /d S e L ('- s >(„ .,,/„ , , , 
7 \d x h(z(s)) - d x g (z(s))J 



o 

We then recall that h(z) satisfies 

\\h{z)\\ 2 a + ||Dfc(z)||,. <q|z|| 2 , 
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which implies 



||M[z]|| 2 3 <Csup(l + t)*B [f](t) <C, ||DM[z]|| 2 | <Csup(l + t)'B [|](«) <C 

4 t>0 4 t>0 

Moreover h (a, b) = f(a, b)d x b + g(a, b) ~~ go(a, b) satisfies 

|Mz)lki + l|DMz)lkr <CN 2 - 



Here, we need to consider separately t E [0,1] and t > 1 when estimating 1 1 PDA/2 [z] 1 1 2 
characteristic function for t > 1, we find 



Writing again Q for the 



\\m 2 [z}\\ 2 3* < c sup gg*^ !■;](«) < c 



t>o 



11(1 



DAT 2 [z]|| 2 5* <C sup (l + t)f <C, 



0<t<l 



Mz]|| 2 -<c S up^4|»J° 0>1 ](t)<c 



t>l 



We finally note that 



and so, using Lemma[T9l we find 



|| 50 (z)|| 2 3 +||D 5o (z)|| 2 5 <C||z|| 



l^aWll a , r <8up^4j;|;; i0 ](t)<c I 



13 1 



t>0 



||DA4N|| 2 , r < ^Pggn 4m 4 o]W ^ c ■ 



f>0 



2 ' 4 ' 2 ' 

1 _ 

2 ' 4 ' 2 ' 



This completes the proof. 



It now only remains to prove the estimates (13.13b and ( 13.14b on the maps TZs u<v \, where we recall that 

TZ u [z,R N }{t) = c + E [h ltU + h 3 , u ]{t) - c-E- 2 [h ltV + h 3 , v ](t) + c 3 E-i[k 2 + h 4 ]{t) , 



with 



Tl v [z, R JV ](i) = c-EolAi,* + h 3<v ]{t) - c+E 2 [h hu + h 3 ^ u ]{t) - cgE^ + h 4 ](t) , 

E a [h](t) =d x ds e d " (t ~ s) T a h{s) and 
Jo 



Here, we will only prove that 



]T \\D a n {UiV} [z,K N }\\ 2d+% _ £ < Ce a J2 ||D a R JV || 2)S+f _ e + C . 



Q=0 



(7.5) 



It is then straightforward to show (13.14b . namely that the maps H{ UjV } are Lipschitz in their second argument; we omit the 
details. 

To prove (17.5b . we first need estimates on hi = {h\ tU , hi tV ), h 2 , h 3 = (h 3)U , h 3>v ) and /14. We note that Uo = (uq, vq) and 
Ui = (ui,vi) satisfy 



|uo||i,o + ||ui||i, + Ipuollx 1 + IpuiHi 1 < C 
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Bup(l+t)*(|uo(±t,t)| + |ui(±t,t)|)+ (1 + <) 2 (|Du (±i,i)l + |Dm(±t,t)|)< C 
for some constant C; see PropositionfT2l We thus find that 

1 

llhilli,!-. + HDhilkj^ + + HD^ 2 |li,3_ e < Ce £ IP^Ikf+f-* . 

a=0 

||h 3 ||i,i + ||Dh 3 ||i,3 + ||Ju||i,a + ||D/l 4 || 2 ,2 < C . 

The proof of (I7.5l l then follows from Proposition[2T| which implies that 

1 l 
£ ||D*E ff [h 1 ]|| 2i|+f e + ||D Q E ff [ft 2 ]|| 2 ,| +f _ e < Ce £ IP^Ikf +* -« . 

1 

£ ||D"E ff [h 3 ]|| 2 ,3 +f . + ||D Q E (T [MII 2 4 +f * < C 

Q=0 

for any a € {—2, — 1, 0, 1, 2} if the estimates in (17.6b are satisfied. 

Proposition 21 Lef e > ami a £ {—2, —1, 0, 1, 2}. 77zen f/zere /Was 

l l 



£ ||D a E ff M| 2! a +f _ 6 < IP Q Mi,H 

a=0 a=0 
1 1 

^||D Q E CT [/ l2 ]|| 2i 3 + f , <c£||D Q fc 2 ||i,H 



a=0 a=0 



Proof. Let itj = E CT [/^]. Taking the Fourier transform, we find 

ft 

Ui{k,t) = ik dse~ k ^-^ +iaks hi(k,s) . 
Jo 

We can restrict ourselves to J2a=o IP^illi.i+f -e = 1 and J2a=o IP^Mki+f = L Then, it follows that 



|D Q Ul || 2 ,3 +f _ e < Csup(l + t)i+f- £ Si[nJ« J] (*) ^ C 
||D^ 2 || 2 , f+r < Csup Mg^i S X [ }+J^ ] (t) < C 



(7.6) 



for a = 0, 1 as claimed. 
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